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(54) Digital cartographic system for geographical information processing. 



(57) In a digital cartographic system for processing geographic information, a map projection of a map is 
converted. The map projection transformation is carried out with a coordinate value in which longitude 
and latitude are set to a center by way of a linear interpolation so as to simplify a calculation required for 
transforming the map projection, and also not depending upon a map projection to be converted. Also, 
after the transformation process operation, a reference line such as longitude/latitude and an 
equiangular line is superposition-displayed/overlapping-displayed in order to assist an analysis, in 
response to a demand. In accordance with the present invention, a map can be shown with correctness 
of a distance and correctness of azimuth by performing various sorts of map-projection transformation, 
the digital cartographic system can support an intuitive idea of a user, and also can improve an 
information representation function of this digital cartographic system. 
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The present invention generally relates to a meth- 
od for displaying graphic data by a designated map 
projection in a digital cartographic system for geo- 
graphical information processing. More specifically, 
the present invention is directed to both of a system 
and a method for varying the displaying map projec- 
tion by transforming a coordinate value of the graphic 
data into a desirable map projection. 

A method for representing shapes, and ups and 
downs of a ground surface of the earth is so-called as 
a map projection (map projection transformations). As 
the necessary conditions of the map projection, there 
are such elements: (J) an actual distance being ana- 
logously represented on a map (equidistance); @ an 
actual area being analogously displayed on a map 
(equivalence); and ® an angle at a ground surface 
being correctly expressed on a map (equiangularity). 
Since a map drawn on paper is such a fact that a 
sphere is projected onto a plane, all of the above- 
described conditions cannot be simultaneously satis- 
fied. To satisfy any of these conditions, a map is rep- 
resented by a proper map projection, depending upon 
its usage. 

The various sorts of map projections are descri- 
bed in detail in, for instance, Japanese publication en- 
titled "Edition and Projection for Map" written by K. Ko- 
saka, published by SANKAIDO publisher in 1982. In 
this publication, the calculations for projecting the co- 
ordinate systems expressed by longitude and latitude 
on a spherical surface onto the coordinate systems on 
a plane, are classified in accordance with shapes of 
projection surfaces and conditions of projection. 
There are employed a plane, a cylinder and a conic 
as the projection planes, whereas the above- 
described three conditions are handled as the projec- 
tion conditions. 

Conventionally, the coordinate transformation in 
case that map data are processed in the computer, is 
handled in view of inputting of a map and calculations 
of a distance and an area. When a map is inputted, it 
is necessary to convert a coordinate value of map 
data which has been written on paper by way of a dif- 
ferent map projection, into a common coordinate sys- 
tem in order to totally manage the data. To manage 
the map data under better matching conditions, it is 
preferable to employ an expression of a coordinate 
system based on the longitude and latitude. If the 
longitude and latitude are utilized, the map data can 
be continuously managed over an entire region of the 
earth. In transactions of Japanese System Control In- 
formation Institute, "A Geographic Information Sys- 
tem REALS for Personal Computer", volume 34, No. 
5, pages 138-146, 1990 by Taniguchi et al., and 
"Summary of Country Numeral Information" issued by 
Japanese Construction Ministry, Geographic map de- 
partment, map supervision division in 1985, data 
management based on the longitude and latitude has 
been performed. In the former publication, there is de- 



scribed such a method for transforming the graphic 
data into the coordinate system by the longitude and 
latitude, which have been inputted by the Universal 
Transverse Mercator coordinate system (UTM coor- 

5 dinate system) employed in a map with a reduction ra- 
tio of 1/25000. Since the UTM map projection corre- 
sponds to an equiangular projection method, and fur- 
thermore a distance can be measured under practi- 
cally acceptable precision in the map with the medium 

10 reduction ratio, there are merits that both the angles 
and distances can be directly calculated from the co- 
ordinate value. As a consequence, the UTM coordin- 
ate system has been employed to calculate the dis- 
tances and areas in the latter publication, and also 

15 there is described such a method for transforming the 
longitude/latitude of the map data into the UTM coor- 
dinate system so as to perform this calculation. 

Among others, there is another publication "Geo- 
graphic Information System with Superior Analyzing 

20 Function: ARC/INFO" written by Imai, PIXEL, No. 54, 
pages 65 to 70, 1987 as a digital cartographic system 
handling a coordinate transformation. This system 
owns the coordinate transformation function as the 
edition function during the input/updating operations 

25 of the graphic data, and then can execute several dif- 
ferent map- projection transformations with respect to 
the map data. 

In case when map data are processed in comput- 
ers, another map-projection transformation is re- 

30 quired in view of output operations other than the 
above-described input operation. It is difficult to di- 
rectly judge both azimuth and a distance from the map 
data which have been totally managed based on the 
longitude and latitude. While the digital cartographic 

35 system is utilized as an information representing 
means to analyze a region, a proper information rep- 
resentation suitable for a desirable analyzing purpose 
is required. To achieve an intuitively understandable 
analysis support, it is required to display map data 

40 with satisfying correctness in azimuth as well as cor- 
rectness in a distance. There are the azimuth map 
projection and the UTM map projection, functioning as 
a method for representing either azimuth, or a dis- 
tance. In accordance with the azimuth map projection 

45 used for a plane projection, azimuth at a center point 
of projection is equal to azimuth at a ground surface, 
and a line segment for connecting this center line and 
an arbitrary point, becomes a minimum path between 
two points. In accordance with the UTM map projec- 

so tion used for the cylindrical-plane projection, an angle 
between a ground surface and a corresponding map 
becomes equal, and also a distance on the map can 
be expressed as being substantially equal to an actual 
distance under practically acceptable map precision 

55 with a medium reduction ratio. However, it is only pos- 
sible to express continuous coordinate system with 
having the interval of longitude within a range of 6 de- 
grees in the UTM map projection. As to an area, it is 
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possible to equally express the areas by way of the 
conical map projection and the cylindrical map projec- 
tion. 

To display map data in various sorts of map pro- 
jections in conformity to objects, coordinate trans- s 
forming operations should be carried out at a high 
speed and also at high precision. Further, when a co- 
ordinate system is converted into such a coordinate 
system as in the UTM coordinate system where dis- 
continuities are present every 6 degrees, the discon- 10 
tinuity should be avoided by a proper way. 

First, there is a problem in the coordinate trans- 
forming speed. Generally speaking, the coordinate 
transforming operation for the map data must be car- 
ried out with respect to a large quantity of graphic 15 
data. In the publication "Edition and Projection of 
Map", the projection formulae from the longitude and 
latitude to the coordinate systems of the various map 
projections are described. Since each of these projec- 
tion formulae involves trigonometric functions and a 20 
logarithm, and therefore requires a large amount of 
calculating steps, there is a problem that a lengthy cal- 
culation time is necessarily required. In particular, 
since the projection formula to the UTM coordinate 
system is expressed by a series expanding formula, 25 
a back projection is not easily performed, but it is not 
suitable to calculate a large amount of calculation ele- 
ments. It should be noted that although the coordinate 
calculating method with employment of the transform- 
ing table between the longitude/latitude and the UTM 30 
coordinate system has been described in the above- 
described publication "Edition and Projection of Map", 
since this transforming table has a high volume in unit 
of 1 minute and also the transforming formula is a bi- 
quadratic polynomial, a total amount of this calcula- 35 
tion is not so reduced. In the publication "Summary of 
Country Numeral Information", in order to simplify the 
transformation from the longitude/latitude into the 
UTM coordinate system, the graphic data are subdi- 
vided into segments and one representative point is 40 
set in the segment, whereby the transforming formula 
is analogous to a local formula. Although the calcula- 
tion of this method may be simplified, there is another 
problem that shifts happen to occur in the graphic data 
on the boundary line over the adjoining segments, re- 45 
suiting in an occurrence of discontinuities. Further- 
more, the coordinate transformation effected in the 
ARC/INFO system is intended to perform the editing 
operation during the input/updating operations of the 
data, but to convert/display the map projections. so 

Subsequently, there is a problem in the disconti- 
nuities of the graphic data. The UTM map projection 
is effective to express an angle and a distance. To 
suppress distortion in the distance, the entire region 
is subdivided into narrow strips called zones at 6 de- 55 
grees in this UTM coordinate system. Therefore, as 
represented in Fig. 8, there are great shifts at a bound- 
ary line of the adjoining zones in accordance with the 



higher latitude, so that the map data are represented 
with having the discontinuities on this boundary line. 
To represent the graphic data extending these zones 
in a continuous form, there is another great problem 
that such a distortion becomes large in accordance 
with the higher latitude. It should be noted that the co- 
ordinate transformations as defined in the above- 
described publications "Summary of Country Numeral 
Information", and "A Geographic Information System 
REALS for Personal Computer", and intended to ef- 
fectively manage the data and also correct the analyz- 
ing amount, but not to represent various information 
to users. As a consequence, these coordinate trans- 
formation methods still contain the above-described 
discontinuities in the graphic representations. 

In addition thereto, it is necessary to readily pro- 
duce reference lines such as longitude, latitude and 
a curve of equal bearing during representation, which 
own such a function to assist the analysis. Features 
of the represented map data may be easily grasped 
by superimposing these reference lines thereon. 

Embodiments of the present invention may pro- 
vide a digital cartographic system for solving the 
above-described problems such as the continuity of 
the graphic representation and the simplification of 
the coordinate transformation calculations, and also a 
method for transforming the coordinate system and 
for displaying the converted coordinate system, and 
furthermore, to a method for displaying a reference 
line having a function to assist an analysis. 

To simplify calculations for transforming a coor- 
dinate system irrelevant to a map projection to be con- 
verted, a transformation table is employed which cor- 
responds to a means for transforming the correspond- 
ing coordinate values, and thus the coordinate trans- 
formation is carried out by the linear interpolation with 
reference to this transformation table in the present in- 
vention. In the transformation table, coordinate valves 
in the various sorts of map projections are related to 
the (equivalent) longitude and latitude and then the 
transformation is carried out via the longitude and lat- 
itude. As a result, it is possible to execute the coordin- 
ate transformation among arbitrary map projections 
contained in the table group. 

Concretely, first, the map data is represented by 
the coordinate system based upon the longitude and 
latitude in order to maintain matching characteristics, 
and the map data is managed in unit of a sufficiently 
small section in order that the coordinate transforma- 
tion by way of the linear interpolation can be achieved. 
As the segment, a region surrounded by longitude 
and latitude separated at a constant interval is util- 
ized. In addition thereto, both longitude/latitude of 4 
vertexes in a segment in unit of respective segments, 
and also the coordinate values at each of the corre- 
sponding map projections are combined as a table 
shown in Fig. 2. If the relationship between the long- 
itude/latitude and the coordinate values at four ver- 
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texes in unit of subdivision management can be found 
out, then a coordinate value at an arbitrary point within 
a segment unit can be calculated from a proportional 
distribution between the vertexes. Also, the reference 
lines such as the longitude, the latitude and the curve 5 
of equal bearing, which can assist an analysis, can be 
calculated from the proportional distribution between 
the vertexes. The transformation is carried out with re- 
spect to the graphic data stored in the database. If the 
coordinate transformation is performed between one 10 
map projection A and the other map projection B, the 
graphic data stored in the database corresponding to 
the representation by the map projection A is referred, 
whereby the coordinate system of the longitude/lati- 
tude thereof is converted into the coordinate system 15 
for the map projection B. It can be prevented deteri- 
oration in positional precision caused by repeatedly 
perform the coordinate transformation from the coor- 
dinate values of the longitude/latitude stored in the da- 
tabase. It should be noted that with respect to a coor- 20 
dinate system for such a map projection that coordin- 
ate values of the respective points are different from 
each other in response to designation, such as an azi- 
muth map projection requiring a base point (namely, 
a center point of projection) and also a conical map 25 
projection depending upon a region designated by a 
user, since a transformation table corresponding to 
the longitude and latitude cannot be previously ex- 
pressed by numeral values, a table forming routine to 
produce a table in response to the coordinate value of 30 
the base point is prepared, and then the transforma- 
tion table is expressed by the numeral values every 
time the base point is given. 

With regard to the UTM coordinate system, a cor- 
responding table between longitude/latitude and UTM 35 
coordinates for a single zone is prepared. The UTM 
coordinate system is periodic. The same coordinate 
values are repeatedly utilized for each zone is shifted 
by 6 degrees along the longitudinal direction. There- 
fore, an overall region of earth can be expressed by 40 
way of a remainder calculation for the longitude. As 
shown in Fig. 9, when a region extending two adjoin- 
ing zones is handled, the discontinuities is avoided by 
property shifting the zone by 6 degrees. 

Since the map pattern is represented in the con- 45 
verted mode in conformity to correctness in a distance 
and correctness in azimuth in accordance with an ob- 
ject by a user, an intuitive idea can be supported. 
Since the linear interpolation is utilized to convert the 
coordinate values of the map pattern, both of the long- so 
itude/latitude and the coordinate values of the respec- 
tive map projections can be mutually and easily con- 
verted therebetween. The transformation operations 
are independently carried out in unit of section by util- 
izing the transformation table corresponding to the 55 
transforming means, so that the calculation process- 
es can be performed in a parallel distribution mode. 
To execute the coordinate transformation by the linear 



interpolation for the vertex data of the section, the 
continuity of the graphic data among the adjoining 
sections can be maintained. Even when the linear ap- 
proximation is carried out by subdividing the calcula- 
tion unit into the small regions, the coordinate values 
can be calculated at high precision. Furthermore, the 
discontinuities occurring between the adjoining zones 
can be avoided by shifting the zone with respect to the 
UTM coordinate system corresponding to the local- 
ized coordinate system. In the drawings: 

Fig. 1 schematically shows a processing opera- 
tion of a map projection transformation embody- 
ing the present invention; 
Fig. 2 schematically represents a relationship 
among coordinate values of vertexes in a region; 
Fig. 3 is a schematic block diagram for represent- 
ing an example of a construction of a computer 
system to embody the present invention; 
Fig. 4 schematically illustrates a construction of a 
program employed in the computer shown in Fig. 
3; 

Fig. 5 is a flow chart for showing an overall proc- 
ess executed in a computer embodying the pres- 
ent invention; 

Fig. 6 schematically shows a coordinate-transfor- 
mation management table; 
Fig. 7 schematically represents both of a map- 
projection management table and a coordinate 
transformation table; 

Fig. 8 is a schematic diagram for showing an UTM 
coordinate system; 

Fig. 9 schematically indicates a method for pre- 
venting discontinuities of the UTM coordinate 
system; 

Fig. 10 schematically shows a relationship be- 
tween a coordinate transforming region and a 
data file; 

Fig. 1 1 is a flow chart (PAD diagram) for explain- 
ing a process to produce a reference line; 
Figs. 12A and 12B are schematic diagrams for 
showing the process for generating the reference 
line; 

Fig. 13 is an illustration for showing a reference- 
line management table; 

Fig. 14 is a flow chart (PAD diagram) for explain- 
ing a process for evaluating precision; 
Fig. 15 is a schematic diagram for explaining a 
process to evaluate precision; and, 
Fig. 16 is an illustration for showing a path man- 
agement table. 

Referring now to Fig. 1 , one embodiment accord- 
ing to the present invention will be summarized. 

In Fig. 1 , map data has been subdivided and man- 
aged in accordance with a coordinate system based 
upon longitude and latitude, and then stored into a 
map graphic database 105. It should be noted in this 
preferred embodiment that transformation/display 
processes are entirely performed via the longitude 
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and latitude coordinate system. In other words, coor- 
dinate data inputted by a user are interpreted as long- 
itude/latitude coordinate data, and coordinate trans- 
formations are executed with respect to the map fig- 
ure data stored in the map figure database 105. As a 
result, it is possible to prevent accumulation of trans- 
formation errors of figure data which have been 
caused during the preceding coordinate transforma- 
tions. 

The transformation process is carried out by a lin- 
ear interpolation manner, so that the transformation 
may be expressed by way of a matrix of 3X3 con- 
stants, and then a highspeed process operation by 
means of a hardware may be achieved. Since the 
transformations may be independently carried out ev- 
ery subdivision management unit to be subject, it is 
easy to execute a parallel distributed processing. It 
should be noted that the transformation display proc- 
ess is performed as follows. Q) The transformation 
conditions are determined every subdivision manage- 
ment unit, whereby the transformations suitable for 
the sections are performed. ® Otherwise, the trans- 
formation condition is determined from a portion to 
which a user pays his specific attention (will be refer- 
red to "a region of interest"), so that the uniform trans- 
formation is performed over the entire region to be 
processed, in accordance with the flexible process 0 
since the transformation process is performed in such 
a manner that the coordinate transformation formula 
is set every section unit, the transformation may be 
executed in uniform precision over the entire region. 
In accordance with the uniform process, since the 
transformation formula is set from the region of inter- 
est and then the resultant transformation formulae are 
applied to the overall region, only one transformation- 
condition setting operation is required, whereby cal- 
culation process workloads for setting the transforma- 
tion conditions may be reduced. 

Both a map projection B to be converted and a re- 
gion to be converted are entered by a user. When a 
base point of transformation is required in this map 
projection B to be converted such as an azimuth map 
projection, this base point is also inputted by the user. 
Both the region and the base point which have been 
entered by the user are stored in a table for managing 
the coordinate transformation. In Fig. 1, a region 1 in- 
dicates a transformation region (by a representation 
map projection A) designated by the user. In case that 
the uniform process is performed, a region of interest 
is furthermore designated. A region 102 surrounded 
by points to P 4 within the region 101 corresponds 
to the region of interest. Although the region of interest 
will be described later, both a coefficient of a coordin- 
ate transformation formula and precision in the trans- 
formation are determined from this region of interest. 

In accordance with the various conditions desig- 
nated by the user, the map projection is converted by 
the system. In case that the region has been entered 



by the user in the coordinate system of the map pro- 
jection A, both the coordinates of the vertexes of the 
region and the coordinate of the base point are back- 
converted from the map projection A into the longitude 

5 and latitude with employment of the transformation in- 
formation utilized when the map is displayed under 
the map projection A. Subsequently, as represented 
in Fig. 2, the coordinate transformation formula is de- 
termined based upon the relationship between the 

10 longitude/latitude coordinate (X, <fr) of the vertex in the 
management unit and the coordinate (X, Y) by the 
map pattern B thereof with respect to each of the sub- 
division management units contained in the region to 
be converted. Then, the graphic data on the region to 

15 be converted are converted and displayed, whereby 
a transformation result 103 is obtained, while sequen- 
tially retrieving the database 105. The setting and 
transforming process of the transformation conditions 
is carried out every single section in the unit of subdi- 

20 vision management. It should be noted that a region 
104 within the display result 103 denotes a display 
range of the region 101 to be converted by the map 
projection B. Furthermore, after the transformation 
process has been completed, superposition dis- 

25 play/overlapping display of a reference line is per- 
formed as an auxiliary operation of an analysis in re- 
sponse to a user's demand. It should be understood 
that a reference line denotes a line such as a longi- 
tude, a latitude, a grid line of an UTM coordinate sys- 

30 tern, and an isometric line around a certain point. Cor- 
rectness in the azimuth of the display result, correct- 
ness in the distance, and graphic distortion when the 
uniform transformation is performed can be intuitively 
represented to users by displaying such a reference 

35 line in the superposition/overlapping display mode. 

Since the transformation processes are indepen- 
dently carried out in unit of section, these transforma- 
tion processes can be executed in a parallel distribu- 
tion manner. In case that a parallel computer having 

40 a hyper cube structure is employed, if each of the sec- 
tions of the management units is allocated to each of 
the calculation processors, the overall process oper- 
ation can be accomplished at a high speed. As a re- 
sult, a large quantity of graphic data can be readily 

45 handled so that the problems can be easily solved. 

In case of such an environment which cannot be 
processed in the parallel distribution processing 
mode, the process operations for reading and trans- 
forming the graphic data need be repeated by the 

so times equal to the quantity of sections in the suitable 
transformation in the unit of management. To avoid 
such a cumbersome repetitive process, the uniform 
process is carried out based upon the region of inter- 
est. It should be noted that since the coordinate trans- 

55 formation is uniformly carried out over the entire re- 
gion, the graphic data outside the region of interest 
does not always succeed to a nature of the map pat- 
tern B. To compensate for this problem, both the rep- 



6 



9 



EP 0 509 839 A2 



10 



resentation of the reference line and the evaluation of 
positional precision in the graphic data are performed 
in accordance with the user's demand after the pre- 
sentation by the map projection B has been accom- 
plished. The precision evaluation is performed by 5 
comparing with each other amounts of shifts among 
the uniformly converted data and the data converted 
every subdivision unit The permissible range in the 
shift amounts (precision in transformation) depends 
upon worse precision in accordance with both of a 10 
size of a region and precision in graphic data which 
have been inputted by a user. 

In Fig. 3, there is shown a construction of a com- 
puter system for carrying and the above-described 
processing operations. A calculation processing ap- 15 
paratus 301 executes ® a coordinate transformation 
calculation, taking account of a transformation map 
projection and a map projection transformation region 
which have been entered via an input apparatus 304 
corresponding to an input device such as a keyboard 20 
and a mouse; @ display data transfer to a graphic dis- 
play apparatus (GT) 303; and also ® controls of over- 
all process operation. In a main storage apparatus 
302, there are stored a program of coordinate trans- 
formation/display, data on a coordinate transforma- 25 
tion table, and figure data. An auxiliary storage appa- 
ratus 310 includes a program region 31 1 , a region 312 
for a coordinate transformation table file, and a region 
3 1 3 for a coordinate transformation table file, and a re- 
g ion 313 for a graphic data file. Each of storage con- 30 
tents is read out from them into the main storage ap- 
paratus 302, if required. Reference numeral 305 indi- 
cates an I/O interface circuit. 

In Fig. 4, there is shown a constructive example 
of the programs which have been prepared in the pro- 35 
gram storage region 31 1 within the auxiliary storage 
apparatus 310. Reference numeral 401 denotes a job 
control program to control executions of the respec- 
tive process programs; reference numeral 402 indi- 
cates a map projection input program for causing a 40 
user to select the map projection B to be converted; 
reference numeral 403 represents a region input pro- 
gram for causing a user to designate a region to be 
converted; and also reference numeral 404 indicates 
a coordinate input program which is used to input a 45 
point and a base point when a region is designated. 
Furthermore, reference numeral 405 represents a 
program for operating a graphic transformation man- 
agement table, by which a management table for re- 
gion information entered by a user is operated. Ref- 50 
ere nee numeral 406 is a program for forming a coor- 
dinate transformation table, which is used for such a ♦ 
map projection that the coordinate transformation ta- 
ble cannot be previously expressed by numeral data. 
As the map projection capable of not expressing the 55 
table as the numeral data, there are, for instance, an 
azimuth map projection requiring a base point, and 
also a conical map projection acceptable for a region 



designated by a user. This table production program 
406 is arranged by subprograms used for the respec- 
tive map projections. Reference numeral 407 indi- 
cates a graphic-data file operation program for oper- 
ating a graphic data file; reference numeral 408 de- 
notes a coordinate transformation execution program 
by way of the linear interpolation, reference numeral 
409 represents a display scaling program for execut- 
ing a scaling operation when transformation graphic 
data is displayed; reference numeral 410 is a preci- 
sion evaluating program for transformation; reference 
numeral 41 1 indicates a reference line producing pro- 
gram used for superposition/overlay displaying the 
reference line; reference numeral 412 is a display 
control program for controlling reswitching of displays 
among the map projections A and B; and reference 
numeral 413 indicates an input/output program used 
to read/write the above-described programs from/into 
the input apparatus 304 and the auxiliary storage ap- 
paratus 310. 

Fig. 5 is a flow chart (PAD diagram) for represent- 
ing one example of an overall process executed in the 
computer. A process flow will now be described in ac- 
cordance with this flow chart. 

At a step 501 , a map display method according to 
the present invention is initiated by a user. As a result, 
the job control program 401 is read out from the pro- 
gram region 31 1 and then is written into the main stor- 
age apparatus 302 by the calculation apparatus 301. 
Since the process operations are similar when the re- 
spective programs are initiated, the explanations 
thereof are omitted. Upon initiation of the job control 
program 401 , the present display map projection and 
the data on the display range are read out from the 
system table into which they have been previously 
registered, and then prepared in the main storage ap- 
paratus 302. Also, the transformation management 
table operation program 405 is initiated, and the Pres- 
ent display map projection data which have been reg- 
istered in the coordinate transformation management 
table. This display map projection data is compared 
with the previously prepared data. If the display map 
projection data is different from the previously pre- 
pared data, the display map projection data is rewrit- 
ten and this data is newly registered as the present 
display map projection data. 

At a step 502, the job control program initiates the 
map projection input program 402 so as to cause a 
user to input the map projection to be converted. As 
one example, there is such a method that the map pro- 
jections to be converted are represented in a menu 
form and a user will select a proper map projection. 

At a step 503, the map projection management ta- 
ble selected by the user is read. The transformation 
management table operating program 405 reads the 
map projection management table from the coordin- 
ate transformation table file region 312 and writes this 
management table into the main storage apparatus 
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302. This map projection management table is utilized 
to confirm whether or not the map projection to be con- 
verted requires the base point, and whether or not the 
coordinate transformation table must be formed, and 
also to refer to the transformation table. 5 

To execute steps 504 to 511 (will be discussed 
later), the job control program 401 initiates the region 
input program 403 and the coordinate input program 
404. 

At the step 504, the region input program 403 10 
causes a user to enter a region to be converted. This 
is done in such a way that a coordinate of a vertex may 
be entered while requesting the user to input the nu- 
meral data from the keyboard, or a region may be des- 
ignated on the display screen by utilizing a mouse. 15 

In response to an instruction issued by a user, co- 
ordinate data may be inputted in a form of a coordin- 
ate value and also values of longitude and latitude by 
the present display map projection. 

At a step 505, the management program confirms 20 
whether the method for performing this transformation 
corresponds to the suitable process or the uniform 
process. In case of the uniform process, a region of 
interest similar to that of the step 504 is inputted by the 
user at a step 506. 25 

At the step 507, in accordance with the map pro- 
jection management table which has been read at the 
step 503, confirmation is made whether or not a base 
point should be inputted. If necessary, a point input 
operation similar to that of the step 504 is carried out so 
by the user at a step 508. 

With the above-described input operations, the 
map projection transformation has been prepared. To 
prepare this map projection transformation, the job 
control program initiates the coordinate transforms- 35 
tion execution program 408. 

Processing steps from 509 to 511 correspond to 
such a process that the coordinate values inputted by 
the user are equal to those used in the present display 
map projection A. First, a judgement is made at the 40 
step 509 whether the user's input operation is caused 
by the map projection A, or the longitude and latitude 
are directly designated by the user. In case of the 
user's input operation caused by the map projection 
A, the process defined at the step 51 0 is executed. At 45 
this step 510, the coordinate transformation table for 
the map projection A and the longitude/latitude is pre- 
pared. Next, at the step 51 1 , the back transformation 
from the map pattern A to the longitude/latitude is car- 
ried out. With respect to the respective vertexes of the so 
region to be converted and the region of interest which 
have been inputted by the user, a back- 
transformation formula is produced based upon the 
coordinate relationship among the subdivision man- 
agement units containing the vertexes thereof, and 55 
then the map pattern A is converted into the longitude 
and latitude coordinate system. Also, the back- 
transformation from the longitude/latitude coordinate 



system to the map pattern A is executed under control 
of the coordinate transformation execution program 
408. 

At a step 512, a coordinate transformation man- 
agement data table is formed which is used in the map 
pattern B to be converted. Although the contents of 
this data table will be described later, the tables indi- 
cative of the information on the base point and the co- 
ordinate relationship as shown in Fig. 2 are prepared 
for the region to be converted and the region of inter- 
est. 

In such a case that the coordinate transformation 
table has not yet been expressed as the numeral val- 
ues by way of the map projection to be converted such 
as the azimuth map projection and the conical map 
projection, a subprogram corresponding to the sub- 
programs included in the coordinate transformation 
table forming program 406 is first initiated. A coordin- 
ate system in the map projection B is calculated under 
control of the program 406 based upon the longi- 
tude/latitude coordinate values with respect to the re- 
spective vertexes of the region to be converted and 
the region of interest 

When the coordinate transformation table has 
been prepared, one vertex of each of these region to 
be converted and region of interest is picked up. It is 
assumed that this picked up vertex is called as a ver- 
tex "a". From the longitude/latitude coordinate system 
of this vertex "a", the management unit containing the 
vertex "a" on the graphic data file region 313 is deter- 
mined. Subsequently, the transformation data about 
the respective vertexes of this management unit are 
read out, and then a corresponding table between the 
longitude/latitude and the coordinate system of the 
map projection B is produced. The coordinate values 
at the vertex "a" by the map projection B are linear-in- 
terpolated under control of the coordinate transforma- 
tion execution program 408, and the linear-interpolat- 
ed coordinate values are registered together with the 
longitude/latitude coordinate values thereof into the 
coordinate transformation management table. The 
above-described process operation will be executed 
with regard to all of the vertexes of the region to be 
converted and of the region of interest. 

At a step 51 3, an enlargement ratio (reduction ra- 
tio) for display purposes is determined. The enlarge- 
ment ratio (reduction ratio) is determined in such a 
manner that the transforming region converted into 
the map projection B can be stored within the display 
area of the screen. 

First, a circumscribed rectangle is obtained with 
respect to the transforming region expressed by the 
coordinate system on the map projection B, and there- 
after the enlargement (reduction) ratio is determined 
in such a manner that the longer side of the rectangle 
can be stored within the display region. This ratio will 
be referred by the display scaling program 409. 

Next, the map projection and representation with 
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regard to the graphic data of the transforming region 
are carried out. The job control program 401 initiates 
the graphic data file operation program 407 and the 
display scaling program 409. 

At a step 514, both of the transformation and rep- 
resentation for the graphic data are carried out The 
graphic data about the subdivision management unit 
group containing the region to be converted are read 
out under control of the graphic data file operation pro- 
gram 407, and are converted under control of the co- 
ordinate transformation execution program 408, and 
thereafter are enlarged/reduced under control of the 
display scaling program 409, whereby the resultant 
graphic data are transferred to the display apparatus 
303. In case that the transformation executing method 
has been performed every subdivision management 
unit, the process operation is carried out every unit of 
section. The graphic data in the unit of section is read 
out, the coordinate transforming formula is conducted 
from the relationship between the coordinate values 
at the vertexes of this section, and thereafter the co- 
ordinate system of the graphic data is converted. If the 
transformation executing method is performed based 
on the region of interest, the coordinate transforming 
formula is conducted from the relationship between 
the coordinate values at the vertexes of the region of 
interest, and then the subdivision units to be convert- 
ed are uniformly converted. 

A method for obtaining a subdivision manage- 
ment unit group in question from a region to be con- 
verted among a graphic data file will be described lat- 
er. 

The above-described operations are the major 
step of the map projection transformation. Thereafter, 
the job control program 401 executes both the preci- 
sion evaluation and the superposition display/over- 
laying display of the reference line such as the longi- 
tude and the latitude in response to user's request. In 
this case, the precision evaluation program 410, the 
reference line generation program 411, and the dis- 
play change control program 412 are initiated in re- 
sponse to the request. 

A step 515 implies that the process operations 
from a step 516 to a step 520 are repeated until a user 
instructs an end of the process operation. 

At a step 51 6, a judgement is made of an instruc- 
tion made by the user, and the subsequent step is de- 
termined. 

At a step 517, the positional precision on the 
graphic data which has been designated by the user 
via the region of interest, is evaluated under control of 
the display scaling program 409, and then the preci- 
sion evaluated result is displayed on the display re- 
gion under control of this program 409. 

At a step 51 8, the reference line is superposition- 
displayed/overlapping-displayed. The longitude, lati- 
tude, and grid line of the UTM coordinate system are 
superposition-displayed/overlapping-displayed, and 



erased at intervals designated by the user. 

At a step 519, the representations between the 
map projection A and the map projection B are re- 
changed. The transformation display executed in this 

5 embodiment is returned to the original map projection 
A, otherwise set to the map projection B. As a conse- 
quence, the display results of both of the map projec- 
tions A and B may be compared with each other, and 
thus the analysis support may be performed more 

10 flexibly. The job control program 401 initiates the dis- 
play control program 412 used for changing the rep- 
resentation, so that the back transformations from the 
map projection B to the longitude and latitude, and 
also from the longitude and latitude to the map projec- 
ts tion A are executed in accordance with the contents 
of the coordinate transformation management data 
table. Otherwise, the back transformation from the 
map projection A to the map projection B is again exe- 
cuted. At this time, since all of parameters required for 

20 the transformations have been prepared, the display 
control program 412 merely controls the executions 
from the coordinate transformation table producing 
program 406 into the display scaling program 409. It 
should be noted that the coordinate transformation ta- 

25 ble producing program 406 executes changing oper- 
ations of the coordinate transformation management 
data between the map projection A and the map pro- 
jection B. When the transformation results are dis- 
played on the separate windows, each of these trans- 

30 formation results is displayed on the separate win- 
dows, and these windows may be changed. More 
specifically, when both of these windows are posi- 
tioned side by side in such a manner that these win- 
dows are not overlapped with each other as perme- 

35 able as possible on the display screen, an easy win- 
dow comparison can be achieved and an easy repre- 
sentation can be realized. 

At a step 520, an end process is executed under 
control of the job control program 401. The display 

40 map projection data to be converted within the coor- 
dinate transformation management table is again reg- 
istered to the present display map projection data by 
the transformation management table operation pro- 
gram 405, and thus the map projection data to be con- 

45 verted is erased. Thereafter, an end instruction for the 
step 516 is given. 

Referring now to Figs. 6 and 7, tables for manag- 
ing information about the map projections and map 
projection transformations will be described. Both of 

50 the information about the map projections and the in- 
formation about the coordinate transformation are in- 
dependently prepared from the transformation proc- 
ess, so that the transformation process may be simply 
executed and the problems may be solved. Fig. 6 and 

55 7 show contents of a coordinate transformation man- 
agement table, a map projection management table 
and a coordinate transformation table. 

Fig. 6 shows the coordinate transformation man- 
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agement table. A coordinate transformation manage- 
ment table 601 corresponds to a table used for estab- 
lishing a relationship among the present display map 
projection, the longitude/latitude management data 
and the map projection to be converted as shown in 
Fig. 1. The coordinate transformation management 
table 601 is constructed of two data tables for combin- 
ing the longitude/latitude data in the graphic data file 
region 313 with the projection data to be converted, 
and also of a pointer table for designating these data 
tables. First, the pointer table is arranged by a pointer 
602 to the present display map projection data table, 
and also a pointer 603 to a display map projection 
data table to be converted. A coordinate transforma- 
tion management data table 610 designated by the re- 
spective pointers, is constructed of an identifier 61 1 of 
a map projection; transformation method information 
612 indicating a preparation state of a transformation 
table and also such a matter whether or not a region 
of interest is designated; a base point coordinate 61 3; 
a vertex coordinate table 614 for a transforming re- 
gion; a vertex coordinate table 615 of a region of in- 
terest; and a pointer 616 for pointing out the coordin- 
ate transformation table. The transformation method 
information 612 indicates discrimination whether the 
transformation is performed based upon the respec- 
tive management unit, or the uniform transformation 
is carried out based on the region of interest; discrim- 
ination whether or not there is a base point; and dis- 
crimination whether or not the coordinate transforma- 
tion table has been prepared with respect to an overall 
region to be converted. The relationships between the 
longitude/latitude coordinate values of the vertexes of 
the region shown in Fig. 2, and the coordinate values 
by the map projection denoted by the identifier 611 
have been registered in the vertex coordinate tables 
614 and 615. When the contents of the coordinate 
transformation management table are rewritten at a 
step 520, merely the pointers are substituted. The 
transformation table pointer 616 is arranged by a poin- 
ter 61 7 used for the longitude/latitude coordinate table 
and a pointer 618 used for the coordinate table of the 
map projection indicated by the identifier 61 1 , and is 
accepted from each of the map projection manage- 
ment tables 701. 

In Fig. 7, there is shown a relationship between 
the map projection management table and the coor- 
dinate transformation table. A map projection man- 
agement table 701 corresponds to such a table for 
managing information related to a map projection 
shown in a map projection identifier 702. This man- 
agement table 701 is constructed of graphic informa- 
tion 703 related to such an information whether or not 
the transformation table has been expressed by nu- 
meral values, and also whether or not the base point 
is required; an identifier 704 of a transformation table 
forming program; and a pointer 705 for pointing out a 
transformation table within the coordinate transforma- 



tion table file region 312. The graphic information 703 
consists of a NUMERAL-LABEL and a BASE-POINT- 
LABEL The NUMERAL-LABEL represents whether a 
transformation table is already expressed as the nu- 

5 meral values or not. The BASE-POINT-LABEL repre- 
sents whether a map projection requires a base point 
or not. For instance, as to the azimuth map projection, 
the NUMERAL-LABEL is set to 'YET". (The NUMER- 
AL-LABEL requires to make the table expressed as 

10 the numeral values.) The BASE-POINT-LABEL is set 
to "NEED". (The azimuth map projection requires a 
base point.) On the other hand, as to the conical map 
projection, the NUMERAL-LABEL is set to "YET and 
the BASE-POINT-LABEL is not set to "NO-NEED". 

15 (The conical map projection does not require a base 
point." If the transformation table is not expressed as 
the numeral values, the transformation table file poin- 
ter 705 is not stored, but the program identifier 704 
has been set instead of this pointer. At a step 513, the 

20 program corresponding to the program identifier 704 
in the coordinate transformation table forming pro- 
gram 406 is initiated. If the transformation table has 
been expressed as the numeral values, the program 
identifier 704 is set to "0", and then a pointer for des- 

25 ignating a file is stored into the transformation table 
file pointer 705. This transformation table file point 
705 is arranged by a pointer 706 for the longitude/lat- 
itude coordinate file and also a pointer 707 for a co- 
. ordinate file of a map projection indicated by the iden- 

30 tifier 704. 

A transformation table file group 71 0 is stored into 
the coordinate transformation table file region 312 in 
accordance with the management unit of the graphic 
data which has been stored into the graphic data file 

35 region 313. The transformation table file group 710 is 
arranged by a longitude/latitude file 71 1 indicative of 
a longitude/latitude coordinate system corresponding 
to each of vertexes in the respective management 
units, and a coordinate value data file 712 provided for 

40 each map projection corresponding to the same ver- 
tex in conformity to the minimum subdivision manage- 
ment unit of the graphic data. In case that the trans- 
formation table is not previously expressed in a nu- 
meral value, these tables 71 1 and 712 of the subdivi- 

45 sion management unit group corresponding to the re- 
gion to be converted are produced with employment 
of the coordinate transformation table fonning pro- 
gram 406. The produced tables are registered at table 
pointers 617 and 618. The coordinate relationship 

50 among the vertexes in the subdivision management 
unit, can be easily conducted from the relationships 
720 among the sections and vertex data, if the subdi- 
vision management unit is determined. This process 
operation is carried out under control of the transfor- 

55 mation management table operation program 405. 

Figs. 8 and 9 are explanatory diagrams in case of 
the UTM coordinate system. 

Fig. 8 is an explanatory diagram of a coordinate 
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transformation table in the UTM coordinate system. 
As a shown in this figure, to suppress distortions in 
distances during projections, the UTM coordinate sys- 
tem is segmented into zones every 6 degrees in a 
longitudinal direction, and thus an overall surface of s 
the earth is represented by 60 zones. Since the UTM 
coordinate system corresponds to such a coordinate 
system every zone, the same coordinate values are 
obtained every time the zone is shifted by 6 degrees 
along the longitudinal direction. As a consequence, it 10 
is not required to prepare transformation tables over 
the entire file of the graphic data file region, and then 
a transformation table between longitude/latitude and 
the UTM coordinate system with respect to a single 
zone may be prepared. The UTM coordinate system is 
can be related to the residue calculation of the longi- 
tude. Furthermore, if an attention is given to symme- 
try, only a transformation table with respect to a por- 
tion of a zone corresponding to 1/4 upper portion 
thereof is prepared. The format of the transformation 20 
table is followed to the files 71 1 and 71 2, whereas the 
relationship between the subdivision unit and the ver- 
tex is followed to the relationship 720. The process 
operations including the remainder calculation of the 
longitude are carried out under control of the transfor- 25 
mation management table operation program 405. 

Fig. 9 indicates a method for preventing a discon- 
tinuity of the UTM coordinate system. A description 
will now be made of such a process executed in that 
a region 901 to be converted extends over two zones 30 
902 and 903. By performing this process, the overall 
surface of the earth is segmented again into another 
UTM zones by shifting zones properly. New zones 
can be set including an old boundary and the region 
901 in one zone. Therefore the region 901 has no 35 
boundary in a new zone, so that the proglems can be 
also solved. If an extension of the region 901 along the 
longitudinal direction is within 6 degrees, a zone such 
as the zone 904 is properly shifted, the discontinuity 
in representation can be avoided. An amount of shift- 40 
ing angle in the longitudinal direction is determined 
based upon a central value of the longitude of the re- 
gion to be converted. A similar process operation can 
be achieved by shifting it by this shift amount while 
performing the residue calculation of the longitude 45 
with reference to the corresponding table. 

Next, a coordinate transformation process oper- 
ation by way of the linear interpolation will now be ex- 
plained. In accordance with the coordinate relation- 
ship shown in Fig. 2 a coordinate value (X i( Yj) in the so 
designated map projection corresponds to a longi- 
tude/ latitude coordinate (X,, <h) of a vertex {X u <h) (note 
that i=1 to 4). When the region to be converted is sur- 
rounded by a longitude and a latitude, since A3 = X2\ " 

= X1, <()2 = <t>1 and <t>3 = <|>4, (X, Y) corresponding to 55 
a point (X, <t>) will be calculated from both proportional 
distributions about the respective "X" and proportional 
distributions about the respective "<t> w as follows: 



X = [(<j>2 - <t>){(X - VI) X2 + {X2 - X)X1} + (<(> - 
<|>1){(X - X1)X3 + (A2 - X) X4}] / {($2 - $1) {XZ 
- X1)> 

Y = [(<|>2 - 4>){(X - VI) Y2 + (X2 - X)Y1) + • 
^){{X - VI) Y3 + (X2 - X) Y4}] / {(<|>2 - W)(X2 
- 

In other words, it may be obtained by the linear trans- 
formation while the above-described formulae are ex- 
panded to acquire the below-mentioned formula 
[X, Y, I] = [X , +, V<M . T 
T: 3X3 transformation matrix (constant matrix) 
Since the transformation matrix T is determined only 
from 4 vertexes of the segment and also the transfor- 
mation formula is expressed by the linear form, the co- 
ordinate transformation calculation may be simplified. 
Furthermore, since the coordinate value after the 
transformation are determined from the proportional 
distributions of the coordinate values of the respective 
vertexes, there is no shift in the graphic data on the 
boundary line of the adjoining segments and thus the 
problems can be solved. Since the subdivision man- 
agement unit corresponds to the segment surrounded 
by the longitude and the latitude, the transformations 
can be executed every segment in the suitable proc- 
ess. To the contrary, in the uniform process, the 
above-described formula may be applied by that eith- 
er the region of interest is inputted as the region sur- 
rounded by the longitude/latitude, or the region of in- 
terest is substituted by the longitude/latitude. As the 
substitution method, it is realized by selecting such a 
region which is inscribed with the region of interest 
and surrounded by the longitude/latitude. 

Fig. 10 represents a relationship between a re- 
gion to be converted and a subdivision management 
unit. A description will now be made of a method for 
obtaining a subdivision management unit group cor- 
responding to a region to be converted from a graphic 
data file. Reference numeral 1001 indicates a data file 
within the graphic data file region 31 3. With respect to 
the respective vertexes of the region to be converted, 
points corresponding thereto on the data file are ob- 
tained from the longitude/latitude coordinates of the 
vertexes. A region 1002 formed by connecting the 
corresponding points of the respective vertexes is as- 
sumed to as a region to be converted on the data file. 
A minimum set 1 003 of a section or segment including 
this region 1002 corresponds to a management unit 
group which should be converted. There is one meth- 
od for managing the minimum set 1 03 as a roster form 
list. The coordinate transformation processes are car- 
ried out with respect to the respective segments. 
Since the transformation processes are independent- 
ly carried out in the unit of segment, these transforma- 
tion processes may be performed in the parallel dis- 
tributed mode. Further, when a transformation table is 
formed, a relation is obtained between the longitudes 
and latitudes at these vertexes of the segments, and 
the coordinate values under the designated map pro- 
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jection. 

Subsequently, a method for generating a refer- 
ence line according to a preferred embodiment of the 
present invention will now be described. This is done 
with respect to a display screen. The representation 
made on this display screen is a region to be convert- 
ed which has been designated by a user. As a result, 
the process to generate the reference line is carried 
out with regard to a graphic data subdivision manage- 
ment group similar to a subject of the coordinate 
transformation process. Here, the kind of reference 
line is designated by the user, the positions of this ref- 
erence tine are calculated every subdivision manage- 
ment unit and displayed. Standard condition values 
are previously set to a display interval and a display 
color of the reference line. A confirmation is made dur- 
ing the setting operation whether or not a setting 
change is requested. If there is no specific input, the 
display of the reference line is carried out in accor- 
dance with the standard conditions. The reference 
line is generated in such a way that a cross point be- 
tween a boundary line (longitude/latitude) of a seg- 
ment and a reference line is calculated, and the cal- 
culated cross points are connected to each other. 
Since the management table of the reference line is 
independently provided from the map-projection 
transformation process and the transformation man- 
agement table, the reference line can be produced ir- 
relevant to the original display map projection and the 
transforming map projection. As a consequence, the 
reference line functioning as an assistance of analysis 
may be freely selected and therefore an analysis sup- 
port may be flexibly achieved. 

Fig. 11 is a flow chart for explaining the process 
to generate the reference line. 

At a step 1 1 01 , a sort of reference line is set. As 
one example, there is such a method to select the sort 
of reference line by a user while the sorts of reference 
lines are displayed in a menu form. 

At a step 1 102, a check is done whether or not a 
management table of the reference line which has 
been set at the previous step 1101 is prepared. If the 
management table has not yet been prepared, a proc- 
ess defined at a step 1103 is executed to prepare 
such a management table. Under control of the sup- 
port-management table operating program 405, the 
reference line management table is read out from the 
coordinate transformation table file region 312 and 
then is written into the major storage apparatus 302. 
Furthermore, under control of the support-manage- 
ment table operating program 405, a map projection 
identifier of the coordinate transformation manage- 
ment table is referred. If either this map projection 
identifier of the coordinate transformation manage- 
ment table is different from another map projection 
identifier of the reference line management table, or 
there is no coordinate transformation table for an 
overall region to be converted, a table is produced in 



accordance with the coordinate transformation table 
producing program 406. A map- projection manage- 
ment table indicated by the map projection identifier 
will be referred, if required. Also when no base point 

5 has been entered in case that the base point is re- 
quired such as in the equiangularity, the base point is 
inputted In accordance with the coordinate input pro- 
gram 404. Upon formation of the table, the resultant 
table is registered in a transformation table pointer of 

10 the reference line management table. If there is the 
same map projection in the coordinate transformation 
management table, it is also registered into a table 
pointer thereof. 

At a step 1 104, an interval of the generation of the 

15 reference line is set. This interval is indicated as Xint, 
Yint and a as shown in Fig. 12. In case of the longi- 
tude/latitude, or the equiangularity, this interval is an 
angle such as every certain degree. In case of the grid 
line of the UTM coordinate system, this interval cor- 

20 responds to a distance interval such as every certain 
kilometers. 

At a step 1105, a display condition of the refer- 
ence line is set. In this preferred embodiment, a dis- 
play color, a width of the reference line, and a sort of 

25 reference line are set. 

A step 1 1 06 implies a present display range. That 
is to say, the process operations defined from the step 
1107 to the step 1115 are repeatedly executed with 
respect to the respective subdivision management 

30 units of the region to be converted which has been 
designated by a user. Since these process operations 
can be separately executed with respect to each of 
the minimum sets 1003 of the segment, the parallel 
process operation may be performed similar to such 

35 a case of the coordinate transformation and thus the 
process operations can be carried out at a high speed. 

At a step 1 107, as shown in Fig. 2, coordinate val- 
ues of the reference amount corresponding to the 
longitude and latitude are combined with each other 

40 with respect to the vertex of the subdivision unit to be 
processed at this time. This method is similar to such 
a method for producing a corresponding table be- 
tween the longitude/latitude and the coordinate of the 
map projection B at the step 512 during the coordinate 

45 transformation. In case of the grid line in the UTM co- 
ordinate system, it may be combined with reference 
to a transformation table by way of the remainder cal- 
culation of the longitude and latitude. 

At a step 1 1 08, a minimum value and also a max- 

50 imum value of a coordinate value for the reference 
amount are obtained from the corresponding table 
produced at the previous step 1107. In case of the 
equiangularity, the coordinate value of the table is 
converted into angular information around the base 

55 point 1208, and then a maximum (minimum) value 
thereof is obtained. A range of the reference amount 
to be processed is restricted by previously obtaining 
the maximum (minimum) value. 
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A step 1109 implies that process operations de- 
fined at a step 1110 and the step 811 are executed 
with respect to each of the reference amounts X (or 
a) contained in between the minimum value and the 
maximum value. At this time, the reference amount is 
expressed as a value obtained by multiplying the in- 
terval value by "K" (symbol "K" being natural num- 
bers). 

At a step 1110, a cross point between the refer- 
ence line and each of the boundary lines (longi- 
tude/latitude) at the subdivision management unit is 
obtained. In case of the longitude and latitude, or the 
grid line of the UTM coordinate, a cross point between 
a straight line X = KXing and the boundary line is ob- 
tained. In case of the equiangularity, a cross point be- 
tween the boundary line and a straight line having an 
angle "Ka" in a designated direction is obtained. At 
this time, a formula of the straight line is expressed is 
follows: 

Y = tan (ka)-(X - Xo) + Yo 
At a step 1 1 1 1 , a line segment to connect the thus 
obtained cross points with each other is displayed. 
The line segment is enlarged/reduced under control of 
the display scaling program 309 and then is transfer- 
red to the display apparatus 203. 

From a step 1 1 12 to a step 1 1 14, a process op- 
eration similar to that defined from the step 1 1 09 to the 
step 1 1 1 1 is carried out with respect to the reference 
amount Y. 

At a step 1115, a subsequent subdivision man- 
agement unit is obtained. 

Fig. 12 schematically shows a process to gener- 
ate a reference line, namely a relationship between 
the reference line and a subdivision management re- 
gion. Fig. 12A represents such a case of longitude 
and latitude, and also a grid line of the UTM coordin- 
ate system, whereas Fig. 12B shows an equiangular- 
ity. Reference numerals 1201 and 1202 indicate 
boundary lines at the subdivision management unit, 
and become longitude and latitude. Dot lines 1203 
and 1204 indicate reference lines indicative of k-Xint; 
reference numerals 1205, 1206 and 1207 denote 
cross points between the respective reference lines 
and the boundary lines; and reference numeral 1208 
denote a base point of an equiangularity. 

Fig. 13 represents a content of a reference line 
management table. Also as to the reference line, 
since the management data is independent from the 
generation process similar to the map projection 
transformation, various sorts of the reference lines 
can be easily handled. The reference line manage- 
ment table 1 301 is constructed of an identifier 1 302 for 
the reference line; a line color 1303 of the reference 
line; a line width 1304; a sort of line 1305; a display 
interval 1306 of the reference line; and an identifier 
1307 of a map projection to handle the reference 
amount; and also a transformation table pointer 1308. 
The transformation table pointer 1308 is similar to the 



transformation table pointer 616 of the coordinate 
transformation management data table. 

Referring now to Figs. 14 and 15, a method for 
evaluating precision, according to one preferred em- 

5 bodiment, will be described in case that a map projec- 
tion has been converted in accordance with a region 
of interest. In this preferred embodiment, a coordinate 
value obtained by the suitable transformation is set as 
a true value, and precision is evaluated by checking 

10 a shift between another coordinate value obtained by 
the uniform transformation and such a coordinate val- 
ue which has been originally obtained by way of the 
suitable transformation. That is to say, the compari- 
son between the first-mentioned coordinate value and 

15 the second mentioned coordinate value, is performed 
at points obtained in the following sequence and posi- 
tioned on the boundary lines in the respective subdi- 
vision management units in the region to be convert- 
ed. 

20 It is assumed that precision is satisfactory until 
the shift between these coordinate values is within an 
allowable range. A region having arbitrary precision 
may be displayed in an emphasis mode by setting the 
allowable error range in accordance with the following 

25 sequence. 

Fig. 14 is a flow chart for explaining a process to 
evaluate precision. From the coordinate transforma- 
tion formula by the uniform transformation, a straight 
line having a constant X and a constant Y is deter- 

30 mined in the coordinate system (X, Y) of the map pro- 
jection B to be converted, a cross point between this 
straight line and a boundary line of a segment, and 
also a shift in a coordinate value at a cross point pos- 
ition is evaluated. Fig. 15 schematically illustrates a 

35 summary of the precision evaluating process. 

At a step 1401 , an allowable range for a positional 
shift amount of graphic data is determined. This allow- 
able range of the shift amount (precision in transfor- 
mation) is determined in accordance with an area of 

40 a region of interest (namely, the region of interest 
owns how many minimum subdivision units) which 
has been inputted by the user, and also precision in 
graphic data. If the coordinate transformation table at 
the region to be converted has not yet been ex- 

45 pressed by numeral values, referring to the coordin- 
ate transformation management data table 610, a ta- 
ble is formed under control of the table forming pro- 
gram 405. At a step 1 402, a straight line (X=X1 , Y=Y1 ) 
having a constant X and a constant Y is determined 

so from a coordinate value (X1, Y1) of a vertex 1501 
positioned at a lower left position of the region of in- 
terest. Subsequently, the precision evaluation is car- 
ried out along these straight lines. In other words, a 
point 1501 is set to a starting point, and then the pre- 

55 cis ion evaluation is executed along a positive direc- 
tion and also a negative direction. At a step 1403, the 
precision evaluation is carried out along the positive 
direction of the straight line X=X1 set at the previous 
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step 1402. The adjacent sections to the straight line 
X=X1 are successively evaluated along the positive 
direction thereof, and it is shown that the process op- 
erations defined from a step 1421 to a step 1427 are 
executed while the allowable precision can be satis- 
fied. It is now assumed that a cross point with the 
boundary line which is separated most remote from 
the region of interest within a range not exceeding the 
allowable precision, is set to be "PX+". 

Reference numeral 1502 indicates "PX+". At a 
step 1404, the precision evaluation is carried out 
along the negative direction of the straight line X=X1 , 
which is similar to that of the step 1403. A cross point 
obtained this precision evaluation is set to be "PX-". 
Reference numeral 1503 indicates "PX-". At a step 
1405, the precision evaluation is performed along the 
positive direction of the straight line Y=Y1, which is 
similar to that of the step 1403. Across point obtained 
atthis evaluation is set to be "PY+". Reference numer- 
al 1504 indicates "PY+". At a step 1406, the precision 
evaluation is carried out along the negative direction 
of the straight line Y=Y1 , which is similar to that of the 
previous step 1403. A cross point obtained during this 
precision evaluation is set to be "PY- n . Reference nu- 
meral 1505 denotes "PY-". At a step 1407, a straight 
line LY+ which passes through the cross point "PX+" 
along the direction of Y=Y1, is obtained, whereas a 
straight line LY- which passes through the cross point 
"PX-" along the direction of Y=Y1 is obtained. At a 
step 1408, a straight line LX+ which passes through 
the cross point PY+ along the direction of X=X1 is ob- 
tained, and also a straight line LX- which passes 
through a cross point w PY- n along the direction of 
X=X1 is obtained. At a step 1410, a region 1509 sur- 
rounded by the four straight lines LX+, LX-, LY+ and 
LY-, namely a region whose vertex is the cross point 
obtained at the step 1409, is represented in the em- 
phasized mode as the region capable of satisfying the 
precision. Reference numeral 1510 indicates a range 
expressed in the emphasized mode on the display 
screen. 

The steps 1420 and 1430 are a flow chart for ex- 
plaining the precision evaluation process operations 
executed from the steps 1403 to 1406. At the step 
1420, the precision evaluation process for the straight 
line X=X1 is performed, whereas at the step 1430, the 
precision evaluation process similar to that of the step 
1420 is performed with respect to the straight line 
Y=Y1. 

A process operation of a step 1421 implies that 
process operations defined from a step 1 422 to a step 
1425 are repeatedly performed unless the shift does 
not exceed the allowable range. At a step 1422, posi- 
tional information on a cross point 1506 between the 
straight line X=X1 and the boundary line is recorded. 
At the beginning, a point 1501 is recorded. After sec- 
ond times, the cross points which have been obtained 
during the previous time are recorded. At a step 1423, 



a subdivision management unit adjacent to the posi- 
tive (negative) direction of the straight line is found. At 
a step 1424, a longitude and a latitude of a cross point 
between the straight line X=X1 and the boundary line 

5 ofthe unit section found at the step 1423 are obtained, 
and also a coordinate system at a cross point by the 
flexible transformation is obtained. Since either the 
longitude, or the latitude is constant on the boundary 
line ofthe segment, both the longitude and latitude of 

10 the cross point can be easily obtained. The coordinate 
value Xr is obtained by way of the linear interpolation 
from the corresponding relationship of the coordin- 
ates at the vertexes of the section. At a step 1425, a 
judgement is made whether or not an absolute value 

15 of (X1-Xr) functioning as the shift amount is calculat- 
ed. If the shift exceeds the allowable range, the cross 
point obtained during the previous time is set to 
PX+(PX-), whereby the process operation is complet- 
ed. Reference numeral 1507 indicates a cross point 

20 exceeding the allowable range. 

Thereafter, a method for carrying out a map-pro- 
jection transformation/representation along a path 
designated by a user, according to a preferred em- 
bodiment, will now be described. This method may be 

25 applied to such a case that display results are differ- 
ent from each other, depending upon a position of a 
base point, as in the azimuth map projection. The 
base points are produced at a constant interval along 
the path inputted by the user, whereby the map pro- 

30 jection transformation is performed and then dis- 
played. If the base point is determined, the map pro- 
jection transformation/display are executable similar 
to that of the previous case. 

Fig. 16 represents contents of a path manage- 

35 ment table. A path management table 1601 is con- 
structed of a coordinate value 1603 of a vertex in a 
path, a quantity of vertex 1602, and also a distance 
1604 up to the next vertex. Reference numeral 1603 
represents a coordinate value of longitude/latitude of 

40 the vertex. Reference numeral 1604 denotes a dis- 
tance up to the subsequently registered vertex, a val- 
ue of which may be calculated either from the UTM co- 
ordinate value, or a distance between two points on a 
sphere. In accordance with a method for generating 

45 a base point of one preferred embodiment, an accu- 
mulated distance among the respective vertexes is 
obtained from the distance 1604 to each vertex, and 
then the points on the path are interpolated every pre- 
determined interval, thereby obtaining a base point. 

so In accordance with the present invention, the 
maps can be displayed in conformity to correctness in 
a distance as well as correctness in an azimuth by 
performing various sorts of map- projection transfor- 
mation, so that these correct maps can intuitively sup- 

55 port user's idea. Also, since the continuities of the 
graphic data about the adjoining sections can be 
maintained by linear- interpolating the vertex data ev- 
ery section in case that the various sorts of map pro- 
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jection transformation, and furthermore, the disconti- 
nuities of the adjoining zones can be avoided by shift- 
ing the zone with respect to the UTM coordinate sys- 
tem corresponding to the local coordinate system, the 
maps can be displayed without causing inconvenient 
feeling to a user. As a result, the information repre- 
senting function of the digital map information system 
can be improved. It is also possible to provide data on 
the basis of analysis. 



Claims 

1. A digital cartographic system for processing geo- 
graphic information comprising: 

means (310) for storing graphic data 
drawn by a first map projection; 

means (304) for designating a second map 
projection; 

means (301) for transforming the graphic 
data drawn by the first map projection into graphic 
data by the second map projection; and 

means (303) for outputting the graphic 
data converted by the second map projection, 
said transforming means (301) being constructed 
of such a means (601, 610, 701, 710) for trans- 
forming the map projection with a coordinate val- 
ue in which longitude and latitude are set as a 
center by a linear interpolation. 

2. A digital cartographic system for processing geo- 
graphic information as claimed in Claim 1 , where- 
in both of means (410) for evaluating precision 
about a result obtained from the map projection 
transformation, and means (410) for displaying a 
specific region in an emphasized mode are em- 
ployed, whereby the precision of the data on the 
transformation result is evaluated, and also such 
a region capable of satisfying predetermined pre- 
cision is displayed in the emphasized mode. 

3. A digital cartographic system for processing geo- 
graphic information as claimed in Claim 1 , where- 
in means (411) for producing a reference line 
such as longitude and latitude is provided, where- 
by when the converted graphic data by the sec- 
ond map projection is displayed on said output 
means, the reference line is superposition-dis- 
played/overlapping-displayed. 

4. A digital cartographic system for processing geo- 
. graphic information, as claimed in Claim 1 , 

wherein means (304) for designating a path by 
coordinate input means is provided, whereby a 
base point of transformation is sequentially gen- 
erated along the path designated by a user, and 
results of the map projection transformation are 
represented on said output means at each of said 



base points. 

5. A digital cartographic system for processing geo- 
graphic information as claimed in Claim 1 , where- 

5 in said transforming means prevents discontinui- 

ties by properly shifting a zone in case that said 
second map projection is converted into such an 
UTM map projection having discontinuous coor- 
dinate values at a boundary of adjoining zones, 

10 and said output means displays a region extend- 
ing two adjoining zones. 

6. A digital cartographic system for processing geo- 
graphic information, comprising: 

is means (312) forsubdividing and managing 

graphic data; and 

means (31 1 ) for performing map projection 
transformation every section in a unit of subdivi- 
tion management, whereby the map projection 

20 transformations are carried out in a parallel dis- 
tributed mode. 

7. A digital cartographic system for processing geo- 
graphic information, as claimed in Claim 6, 

25 wherein the transformation process is performed 
on the graphic data with a database correspond- 
ing to graphic data about a designated region of 
interest. 

30 8. A digital cartographic system for processing geo- 
graphic information, as claimed in Claim 6, 
wherein means for subdividing and managing 
said map projection is separately provided with 
means for executing said map projection transfor- 

35 mation. 

9. A method for displaying a digital map, comprising 
the steps of: 

designating (502) a second map projection 
40 in order to convert graphic data drawn by a first 

map projection into graphic data drawn by the 
second map projection; 

transforming (509-511) for transforming 
the graphic data drawn by said first map projec- 
ts tion into the graphic data drawn by said second 
map projection; and 

outputting/displaying (514) map data by 
said converted second map projection, wherein 
said transformation is performed with a coordin- 
so ate value in which longitude and latitude are set 
as a center by a linear interpolation (512, 408). 

10. A method for displaying a digital map as claimed 
in Claim 9, wherein precision of data on the trans- 

55 formation results is evaluated (517, 1403, 1404, 
1405, 1406), and a region capable of satisfying 
predetermined precision is displayed in an em- 
phasized mode (517, 1410). 
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1 1. A method for displaying a digital map as claimed 
in Claim 9, wherein when the graphic data drawn by 
said converted second map projection is outputted 
and displayed, a reference line is superposition-dis- 
played/overlapping-displayed(518, 1101-1115). 5 

12. A method for displaying a digital map as claimed 
in Claim 9, wherein a base point of transformation 
is sequentially generated along the path desig- 
nated by a user (507, 508) and results of are rep- 10 
resented on said output means at each of said 
based points. 

13. A method for displaying a digital map as claimed 

in Claim 9, wherein said transforming means pre- 15 
vents discontinuities by properly shifting a zone in 
case that said second map projection is convert- 
ed into such an UTM map projection having dis- 
continuous coordinate values at a boundary of 
adjoining zones, and said output means displays 20 
a region extending two adjoining zones. 
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